MADPH-00-1160 
Fermilab-PUB 00/049-T 
AMES-HET 00-01 

March 2000 



Neutrino Oscillations at an Entry-Level 
Neutrino Factory and Beyond 

O", 

O ! V. Barger 1 , S. Geer 2 , R. Raja 2 , and K. Whisnant 2 

O 

1 Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
2 Fermi National Accelerator Laboratory, P.O. Box 500, Batavia, IL 60510, USA 
^Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA 

OS 

^j- ! Abstract 
oo 

£^ ' We consider the parameters of an entry-level neutrino factory designed to 

■ make the first observation of v e — > v„ oscillations, measure the corresponding 

amplitude sin 2 2#i3 , and determine the sign of the atmospheric-scale 5m 2 2 via 
matter effects. A 50 kt detector, a stored muon energy £7„ > 20 GeV and 10 19 
muon decays would enable these goals to be met provided sin 2 29 13 > 0.01. 
The determination of the sign of 5m 2 2 also requires a baseline L > 2000 km. 
An upgraded neutrino factory with O(10 20 ) decays would enable the first 
observation of v e — ► v T oscillations. With O(10 21 ) decays the effects of a large 
^ , CP-phase could be measured in the case of the large angle matter oscillation 

solution to the solar neutrino anomaly. Our analysis includes a family of 
three-neutrino models that can account for the atmospheric and solar neutrino 
oscillation indications. 



I. INTRODUCTION 



The observation of a deficit of muon-neutrinos in atmospheric neutrino experiments [JT^ 
has paved the way for a new generation of experiments studying neutrino masses and mixing. 
The neutrino sector offers exceptional opportunities for studying some of the most funda- 
mental issues in particle physics, such as the origin of masses and CP violation. A major 
advantage over the quark sector is that neutrino phenomena are free of the complications of 
strong interactions. A comprehensive knowledge of the neutrino mixing matrix may yield 
clues to the old puzzle of why there is more than one lepton family. 

The SuperKamiokande (SuperK) collaboration |]T| has found that the atmospheric neu- 
trino deficit is dependent on L/E u , with greater suppression of the flux with increasing 
L/E u . Moreover, the electron- neutrino rate is L/E v independent and consistent with the 
calculated v e flux. The natural interpretation of the atmospheric data is in terms of — > v T 
oscillations, with maximal or near-maximal mixing and a neutrino mass-squared difference 
^ m atm — 3 x 10~ 3 eV 2 . This is supported by SuperK measurements of the zenith angle 
dependence, which due to matter effects is different for — ► v T and — ► u s , and by n° 
production in neutral current events, which taken together rule out — > v s at 99% C.L [[|. 

There are other indications of neutrino oscillation phenomena at 6m 2 values distinct 
from the Sm 2 tm scale. A long-standing puzzle is the observed deficits of solar neutrinos 0] 
compared to the flux predictions of the Standard Solar Model . There are four regions || 
of oscillation parameters Sm 2 olar , sin 2 26 l solar that can accommodate the present solar data. 
Three of the solutions involve resonance enhancements f7]-pl due to the coherent scattering 
of v e from the dense solar medium. The Large Angle Matter (LAM) solution has Sm 2 olaj: ~ 
3 x l(T 5 eV 2 , the Small Angle Matter (SAM) solution has 5m 2 soUr ~ 5 x 10" 6 eV 2 , and 
the Long Oscillation Wavelength (LOW) solution has Sm 2 olar ~ 10~ 7 eV 2 and large mixing. 
Vacuum Oscillation (VO) solutions [|T2"|,|T3f have 5m 2 ~ !CT 10 eV 2 and large mixing. The 



Sudbury Neutrino Observatory (SNO) experiment [14] in progress may be able to exclude 
some of these solar solutions. 

In addition there is some evidence for — > v e and — > v e oscillations from an ac- 
celerator experiment (LSND) at Los Alamos JIB]. The observed event rates correspond to 
sin 2 2# LSND ~ 1CT 2 with ^^lsnd ~ leV 2 ; a sizeable range of Sm 2 values above 0.1 eV 2 is 
actually allowed. The mini-BooNE experiment at Fermilab | ]L6|| , scheduled to start collect- 
ing data in December 2001 with first results expected by 2003, will determine whether the 
LSND effect is real. 

In the next phase of neutrino oscillation studies long-baseline experiments are expected 
to confirm the v ^ — > disappearance oscillations at the Sm 2 tm scale. The K2K experiment 
[I7| from KEK to SuperK, with a baseline of L ~ 250 km and a mean neutrino energy 
of (E v ) ~ 1.4 GeV is underway. The MINOS experiment from Fermilab to Soudan 



with a longer baseline L ~ 730 km and higher mean energies (E v ) = 3, 6 or 12 GeV, 



is under construction and the ICANOE JT9| and OPERA [^OJ experiments, with (E u ) 



17 GeV and baselines L ~ 730 km from CERN to Gran Sasso, have been proposed. The 
various experiments with dominant and beams will securely establish the oscillation 
phenomena and may measure Sm 2 tm to a precision of order 10% Experiments with 

higher mean neutrino energies should be able to observe r production. 

Further exploration of the neutrino mixing and mass-squared parameter space will re- 
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quire higher intensity neutrino beams and 



beams along with z/ M , 



To provide 



these neutrino beams, muon storage rings have been proposed p2| , P3| in which the muons 
decay in a long straight neutrino beam-forming section, and the muons are produced by 
a muon-collider type muon source [24j|. These "neutrino factories" are now under serious 
consideration p2|-pl|1. The resulting neutrino beams would be sufficiently intense to produce 
thousands of oscillation events in a reasonably sized detector (10-50 kt) at distances up to 
the Earth's diameter [ 2~3| . Some initial studies have been made of the physics capabilities of 
such machines p5|-pT| as a function of the stored muon energy E^, baseline L, and intensity 
/. The focus of the present paper is how to best choose E^, L and / to maximize the physics 
output at an entry-level neutrino factory (hereafter referred to as ENuF) and beyond. 

The present work expands on previous studies in several ways. First, we consider the 
minimal E^ and I needed to accomplish our physics goals. Second, we consider the conse- 
quences of a number of model scenarios that can accommodate the atmospheric and solar 
oscillation indications. Third, we investigate possibilities for measuring CP- violating phases 
25| , |30|j3^ - p5T , p7| . Finally, we investigate possibilities for observing v e — > u T oscillations. 



II. THEORETICAL OVERVIEW 



A. Oscillation Formalism 



The neutrino flavor eigenstates v a are related to the mass eigenstates Uj in vacuum by 
unitary matrix U, 



(1) 



The effect of matter on v e beams §-|rT[|T],|8|^| has important consequences for long- 
baseline experiments. The propagation through matter is described by the evolution equa- 
tion 



. d\v a ) 
dx 



J2J2t^t [SmjjUajUpj + A5 ae 5p e ] \v p ) 



(2) 



where x = ct and A/(2E U ) is the amplitude for coherent scattering of v e on electrons, and 

A = 2V2 G F Y eP E u = 1.52 x 10- 4 eV 2 y e p(g/cm 3 )£(GeV) . (3) 

Here Y e (x) is the electron fraction and p(x) is the matter density. The neutrino oscillation 
probabilities are then P(u a — > vp) = \ (vp{x = L) \ v a {x = 0)) | 2 . 

We solve this evolution equation numerically taking into account the x-dependence of 
the density using the Preliminary Reference Earth model Hlj . We have found that for L less 



than about 3000 km (in which the entire neutrino path is in the upper mantle and the density 
is approximately constant), the results of the exact propagation and those obtained assuming 
constant density agree to within a few percent. However, for larger L (where the neutrino 
path partially traverses the lower mantle and the density is no longer nearly constant) the 
assumption of constant density is no longer valid. For example, for L = 7332 km (the 
Fermilab to Gran Sasso distance), event rate predictions assuming constant density can 
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be wrong by as much as 40%. We always use the numerical solution of Eq. (fj) in our 
calculations. 

For three neutrinos (with a = e, fi,T and j = 1, 2, 3) the Maki-Nakagawa-Sakata (MNS) 
| 4~2| mixing matrix will be parameterized by 

C13C12 C13S12 s 13 e~ tS 

U = j — C23S12 — S 13 S23Ci2e i<S C23C12 — Sl3 s 23 s 12e 4<S C13S23 j . (4) 
S23 s 12 ~~ •5i3C23Ci2e* 5 — S23 c 12 — Si^sS^e 15 C13C23 

where Cjk = cosOjk, Sjk = sin^, and 5 is the non-conserving phase. Two additional 
diagonal phases are present in U for Majorana neutrinos, but these do not affect oscillation 
probabilities. 



B. Measuring #13 and the Sign of 8m\ 2 



The oscillation channels v P 



and v e 



can be explored for the first time at a 



neutrino factory. In addition to a first observation of these transitions, the mixing angle #13 
can be measured, the sign of Sm 2 tm can be determined from matter effects, and the CP phase 
5 could be measured or bounded. With this information, models of oscillation phenomena 
can be tested and discriminated. 

The charged current (CC) interactions resulting from v e —>■ and v e — > oscillations 
produce "wrong-sign" muons (muons of opposite charge from the neutrinos in the beam). In 
the leading oscillation approximation the probability for v e — > in 3-neutrino oscillations 
through matter of constant density is |8|,j38|,j3^ 



P{ye -> v,j) = s 2 23 sin 2 sin 2 A™ , 



(5) 



where 



sin 2 29™ 



sin 2 26 
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G4 -cos 26*) + sin 2 20 



(6) 
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and 



1 32 



1.275m 2 2 (eV 2 )L(km) 
E v (GeV) 



\ 



A 



5m 2 



- cos 26* + sin 2 26 



13 • 



(7) 



32 



Here A is the matter amplitude of Eq. @. Thus even with matter effects the v e — > 
probability is approximately proportional to sin 2 26*13. The experimental sensitivity of the 
v e - * measurements therefore changes almost linearly with sin 2 26*13. 

For v e 
and A ~ 5m 2 



oscillations, the sign of A is reversed in Eqs. (|6|) and (0). For sin 26* <C 1 

:12 > °> Pfa 



f M ) is enhanced and P(u e 



is suppressed by matter effects; 



the converse is the case for —A 



5m 2 2 < 0. Thus a comparison of the v e 



and v e 



CC rates gives information on the sign of Sm 



32- 
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C. v e — > v T oscillations 

The availability of u e and u e beams from a neutrino factory would allow a search for 
u e — > z/ T oscillations. In the leading oscillation approximation the probability for u e — > u T 
oscillations through matter of constant density is |8|, j38| , j39[l 

P(y e _> zy r ) = sin 2 2 &™ sin 2 A™ , (8) 

Thus the u e — ► z/ T probability in matter is also approximately proportional to sin 2 2# 13 . 



D. CP Violation 

In vacuum, CP violation in the lepton sector can be explored by comparing oscillation 
probabilities involving neutrinos with the corresponding probabilities for oscillations involv- 
ing antineutrinos. For three-neutrino oscillations in a vacuum, the probability difference 
is 

P(v a -> up) - P{y a -> up) = ^4 J(sin 2A 32 + sin 2A 21 + sin 2A 13 ) , (9) 



where A jk = 1.275m 2 jk (eV )L(km)/E v (GeV) and J is the CP-violating invariant p|,p 



which can be defined as J = Im{U e 2U* 3 U* 2 U fJ-S }. The minus (plus) sign in Eq. (^j) is used 
when a and f3 are in cyclic (anticyclic) order, where cyclic order is defined as e//r. For the 
mixing matrix in Eq. @, 

J = — sin 2^12 sin 2#i 3 sin 2# 2 3 cos 6>i 3 sin 5 . (10) 
8 

Thus even for 5 = ±90°, J will be small when the #12 , #13 mixing angles are small. 

For |5j7i| 2 | ^ l^ m 2il> the CP- violating probability difference for u e — > u^ is given ap- 
proximately by 

\ r>, \ at - /2.54(5m 2 1 (eV 2 )P(km)\ 
P(u e -+ - P(u e -+ „„) - -4Jsm ^ £ ( l GeV j } 1 J , (ID 



It is evident from Eq. ( pT|) that CP violation is only appreciable in vacuum when the sub- 
leading oscillations (in this case oscillations due to dm 2 ,^ begin to develop |j33| . The same 
qualitative results are expected when neutrinos propagate through matter, although the 
oscillation probabilities are changed. 



III. FAMILY OF SCENARIOS 

With three neutrinos, there are only two distinct 5m 2 values. The evidence for atmo- 
spheric, solar and accelerator neutrino oscillations at three different Sm 2 scales cannot be 
simultaneously accommodated in a three-neutrino framework. Here we set the accelerator 
evidence aside and use three-neutrino oscillations to explain atmospheric and solar data; the 
oscillation scale of the accelerator data is more relevant to short-baseline experiments. 
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TABLE I. Representative neutrino oscillation scenarios. 





1 &<m ^ 1 

\om 32 \ 


1 A'JT) ^ 1 

l om 21 1 


sin 2 

Sill -^"23 


sin 2 9fl,„ 
sill ^"\2 


sin 2 9fl, „ 

bill Z(7^3 


J 1 sin 6 


1) LAM 


(atmos) 
3.5 x KT 3 


(solar) 
5 x 10" 5 


(atmos) 


(solar) 
0.8 


0.04 


0.022 


2) SAM 


3.5 x 10~ 3 


1(T 5 


1 


0.01 


0.04 


0.0025 


3) LOW 


3.5 x 1(T 3 


io- 7 


1 


0.9 


0.04 


0.024 


4) BIMAX 


3.5 x 10~ 3 


5 x 10" 5 


1 


1 









A family of representative scenarios in Table |, denned in the ongoing Fermilab long- 



baseline workshop study f45| , will be adopted for our subsequent analysis; the central value 
l^ m 32l = 3-5 x lCT 3 eV 2 is based on the published SuperK data Q. With further data 
accumulation, a slightly lower central value of 2.8 x 10~ 3 eV 2 is indicated ||, with a 90% C.L. 
range of 2 — 5 x 10 -3 eV 2 . The forms of the mixing matrix U in these scenarios are 



[/(LAM) 



[/(SAM) 



[/(LOW) 



[/(BIMAX) 



/ 0.846 
-0.372 - 0.060e w 

V 0.372 - 0.060e il5 

/ 0.994 

-0.035 - 0.071e i(5 
\ 0.035 - 0.071e i<5 

/ 0.807 
-0.413 - 0.058e w 

V 0.413 - 0.058e i<5 

(3* ^ o \ 



0.523 
0.602 - 0.037e i<5 
-0.602 - 0.037e i5 

0.050 
0.706 - 0.004e il5 
-0.706 - 0.004e i5 

0.582 
0.574 - 0.042e i<5 
-0.574 - 0.042e i5 




i 

_1 

2 

1 

2 



1 

V2 
I 

2 

2 V2 I 



1 

V2 



(12) 



(13) 



(14) 



(15) 



Since sin 2#i3 is not well-known, 
sin 2 26 13 = 0.004, for which 



sometimes we will consider the LAM solution with 



[/(LAM') 



0.850 0.525 
-0.372 - 0M9e iS 0.602 - 0.012e i5 
0.372 - 0.019e itf -0.602 - 0.012e w 



0M2e- iS 
0.707 
0.707 



(16) 



Scenarios 1-3 in Table | represent three-neutrino oscillation explanations of the atmo- 
spheric and solar deficits, with the LAM, SAM, and LOW solar options, respectively. We 
do not address the VO solar solution separately, since the sub-leading Sm 2 effects will not 
be significant and the LOW and VO scenarios will be indistinguishable. Scenario 4 with 
bimaximal atmospheric and LAM mixing [fPJ is interesting because the leading oscillation 
decouples in the v e — > channel (U e3 = 0) and the sub-leading oscillations will therefore 
be more visible. However, in this scenario there are no matter effects on u e propagation and 
the sign of 5m 2 cannot be measured; also CP will be conserved. 

The size of the Jarlskog invariant (modulo sin S) is also shown in Table Q; it is largest for 
the LAM and LOW scenarios (which have only one small angle), smaller for SAM (which 
has two small angles), and vanishes for the BIMAX scenario (in which one angle is zero). 
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Since the observability of CP violation depends on both J and the size of the sub-leading 
oscillation scale 8m 21 (see Eq. (0)), one expects appreciable CP violation only in the LAM 
scenario. 

The neutrino mass ordering can in principle be determined by the effects of matter on 
the leading electron neutrino oscillation ||29|| . To illustrate this, Fig. [1] shows the two possible 
three-neutrino mass patterns. Figure [l]a, with one large mass m.3, has atmospheric neutrino 
oscillations with 8m\ 2 > 0. Figure [TJb, with two large masses m 2 ,m 1 , has 8m\ 2 < 0. For 
scenario 4, with no v e participation in the leading oscillations, there are no matter effects at 
the 8m\ 2 scale to determine the sign of 8m\ 2 . 



IV. AN ENTRY-LEVEL NEUTRINO FACTORY 

Neutrino factories require the development of new accelerator sub-systems which are 
technically challenging. The R&D required for a full-intensity muon source might take 
many years. It is reasonable to consider a strategy in which the R&D needed for the first 
neutrino factory is minimized by building, as a first step, a muon source that provides just 
enough muon decays per year to make contact with the interesting physics. If we also wish to 
minimize the cost of an entry-level facility, we must minimize the muon acceleration system 
and hence the energy of the muons decaying within the storage ring. In this section we 
consider, within the framework of the scenarios listed in Table I, the minimum muon energy 
and beam intensity needed at an ENuF. 

We begin by defining our entry-level physics goal. We take this to be the first observation 
of v e — > Up oscillations at the 10 event per year level. The signal will be the appearance of 
CC interactions which are tagged by a wrong-sign muon. To identify signal events we must 
be able to identify muons and measure their charge in the presence of the accompanying 
hadronic shower from the remnants of the target nucleon. Muons can only be cleanly identi- 
fied and measured if their energy exceeds a threshold E min , which in practice is expected |47 



to be a few GeV. This places an effective lower bound on the acceptable energy of the muon 
storage ring. To illustrate this, consider the v e — > signal in a detector that is 2800 km 
downstream of a 20 GeV neutrino factory. The predicted measured energy distributions 
for CC events tagged by wrong-sign muons are shown for the LAM scenario (Table I) in 
Fig. as a function of E min . As E min increases the signal efficiency decreases. For example, 
with E min = 2 (4) [6] GeV the resulting signal loss is 18% (36%) [55%]. In addition, as 
E m[n increases the measured signal distributions become increasingly biased towards higher 
energies, and the information on the oscillations encoded in the energy distribution is lost. 
We conclude from Fig. |2| that with a 20 GeV storage ring we can probably tolerate an E min 
of a few GeV, but would not want to decrease the storage ring energy below 20 GeV. Hence, 
we will adopt 20 GeV as the minimum storage ring energy for an ENuF. In the following, 
our calculations include a muon threshold E m i n = 4 GeV, and for simplicity we assume the 
detection efficiency is for signal events with E^ < E m i n and 1 for E^ > E m i n . 

We next consider the muon beam intensity required to meet our entry-level physics goal 
if the storage ring energy is 20 GeV. To minimize the required beam intensity we must 
maximize the detector mass (M). Recently 50 kt has been considered |TJ] as a plausible 
although ambitious M. We therefore choose M = 50 kt. In addition to the neutrino 
factory beam properties and the detector mass, the signal event rate will depend upon the 
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baseline and the oscillation parameters. Since, to a good approximation, the signal rate 
is proportional to sin 2 2# 13 , it is useful to define the sin 2 2# 13 "reach" for an experiment as 
the value of sin 2 2# 13 for which our physics goal (in this case the observation of 10 signal 
events per year) will be met. Setting the number of useful muon decays per year to 10 19 , 
the sin 2 2#i 3 reach is shown in Fig. | as a function of baseline and 5m\ 2 with the other 
oscillation parameters corresponding to the LAM scenario. The calculational methods are 
described in Ref . PB| , |2T| . The sin 2 2#i 3 reach degrades slowly as L increases, and improves 
with increasing |5m| 2 |, varying by about a factor of 5 over the 5m 2 2 range currently favored 
by the SuperK results. Note that, for 5m\ 2 in the center of the SuperK range, our entry level 
goal would be met with a 20 GeV storage ring and 10 19 decays per year provided sin 2 29 13 
exceeds approximately 0.01, which is more than an order of magnitude below the currently 
excluded region. 

We now consider how the muon beam intensity required to meet our entry-level physics 
goal varies with the storage ring energy. We will choose a baseline of 2800 km, motivated 
by a consideration of the physics sensitivity of an upgraded ENuF (see next section). The 
number of muon decays required to meet our goal is shown in Fig. |] versus the muon storage 
ring energy for the LAM, SAM, LOW, and BIMAX oscillation scenarios. Note that: 

(i) The energy dependent intensities needed for the SAM and LOW scenarios are indistin- 

guishable. 

(ii) Due to the contributions from sub-leading oscillations, the intensity needed for the 
LAM scenario is slightly less than needed for the SAM and LOW scenarios. 

(iii) With a 20 GeV storage ring 2 x 10 18 muon decays per year would meet our entry level 
physics goals for the LAM, SAM, and LOW scenarios. The dependence of the required 
muon intensity I on the storage ring energy is approximately given by / oc E~ ls . 

(iv) For the BIMAX scenario in which sin 2 2# 13 = only the sub-leading Sm 2 scale con- 
tributes to the signal. With a 20 GeV storage ring a few x 10 20 muon decays per year 
would be needed to observe u e — > oscillations. Although this scenario would be bad 
news for a low-intensity neutrino factory, oscillations driven by the sub-leading 5m 2 
scale might be studied with a higher intensity muon source. 

It is straightforward to use the curves in Fig. |] to infer the intensity required to meet our 
entry level goals for LAM, SAM, and LOW-type scenarios with values of sin 2 2# 13 other than 
0.04. For example, if sin 2 2# 13 = 0.01 (a factor of 20 below the currently excluded value, 
and a factor of 4 below the value used for the curves in Fig. ^) we must multiply the beam 
intensity indicated in Fig. f| by a factor of 4 to achieve our entry-level goal. Provided a 50 kt 
detector with good signal efficiency is practical, within the framework of LAM, SAM and 
LOW-type scenarios, a 20 GeV storage ring in which there are 10 19 muon decays per year in 
the beam forming straight section would enable the first observation of v e — > oscillations 
provided sin 2 2# 13 > 0.01. Note that if sin 2 2# 13 > 0.01 the next generation of long baseline 
accelerator experiments (e.g. MINOS) are expected make a first observation of — > v e 
oscillations. These experiments would not be able to measure matter effects and determine 
the sign of Sm 2 2 . Furthermore, if sin 2 2# 13 <C 0.01 then with several years of running the 
entry level neutrino factory data could be used to place a limit on this parameter about an 
order of magnitude below the MINOS limit. 



S 



Consider next the prospects for exploiting u e — > Vp measurements to determine the 



sign of 5m^ 2 . In ref. |29[ we have shown in the LAM scenario that the sign of 5m^ 2 can 



be determined by comparing the wrong-sign muon rates and/or the associated CC event 
energy distributions when respectively positive and negative muons are stored in the ring. 
The most sensitive technique to discriminate 5m\ 2 > from 5m\ 2 < would be to take data 
when there were alternately positive and negative muons stored in the neutrino factory, and 
measure the resulting wrong-sign muon event energy distributions together with the Up — > Vp 
and Up — > Up event energy distributions. The four distributions can then be simultaneously 
fit with the oscillation parameters 5m\ 2 (including its sign), sin 2 2#i3 , and sin 2 2(9 23 left as 



free parameters ||48|| . In the following we take a simpler approach to demonstrate that, 
provided L is large enough, a neutrino factory that permitted the observation of 10 v e — > Vp 
events per year would also enable the sign of 8m\ 2 to be determined. 
We begin by defining the ratio: 

^ ~ We - ^ ( ^ 

R e p is just the ratio of wrong-sign muon rates when respectively negative and positive muons 
are stored in the neutrino factory. Figure [| shows R e p as a function of L for Ep = 20 GeV 
and 5m\ 2 = ±3.5 x 10~ 3 eV 2 /c 4 . Note that when L > 2000 km the ratio R e p for positive 8m\ 2 
is more than a factor of 5 greater than the value for negative 5mJ 2 . As an example, consider 
a 50 kt detector 2800 km downstream of a 20 GeV storage ring in which there are 10 19 
muon decays per year in the beam forming straight section. Suppose that sin 2 26> 13 = 0.01, 
and assume we know that |5m 2 2 | ~ 3.5 x 10~ 3 eV 2 /c 4 from Vp — > Vp measurements, for 
example. If we store positive muons in the neutrino factory after one year we would expect 
to observe 11 wrong-sign muon events if 5m\ 2 > but only 2 events if 5m 2 2 < 0. To reduce 
the uncertainties due to the lack of precise knowledge of the other oscillation parameters, we 
can then take data with negative muons stored. We would then expect to observe less than 
2 wrong-sign muon events per year if 5m\ 2 > 0, but 6 events per year if <5m 2 2 < 0. Clearly 
with these statistics and several years of data taking the sign of 8m\ 2 could be established. 
From this example we conclude that with a few years of data taking a neutrino factory 
that enabled the observation of 10 Vp — > Vp events per year would also enable the sign of 
<5m 2 2 to be determined provided the baseline was sufficiently long (L > 2000 km) so that 
the prediction for R e p changes by a large factor (> 5) when the assumed sign of 5m\ 2 is 
changed. Hence, provided sin 2 2^13 > 0.01, our entry-level neutrino factory would make the 
first observation of u e — > Up oscillations, measure sin 2 2#i3, and determine the pattern of 
neutrino masses. 



V. ONE STEP BEYOND AN ENTRY LEVEL NEUTRINO FACTORY 

An entry- level neutrino factory is attractive if there is a beam intensity and/or energy 
upgrade path that enables a more comprehensive physics program beyond the initial ob- 
servation of v e — > Vp oscillations and determination of the sign of 5m\ 2 . In this section 
we consider the energy and/or intensity upgrades needed to achieve a reasonable upgrade 
physics goal, which we take to be the first observation of u e — > v T oscillations. 
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In ref. |28| we have shown that in a LAM-like scenario the v e — > v T oscillation event 
rates in a multi-kt detector are expected to be significant if is 20 GeV or greater. In 
the following we consider the neutrino factory beam energy and intensity needed to make 
a first observation of v e — > v T oscillations at the 10 event level in one year of running with 
a fully efficient detector (or several years with a realistic detector). Our detector must 
be able to measure r appearance and, to separate the signal from — > u T oscillation 
backgrounds, measure the sign of the charge of the r. Hybrid emulsion detectors in an 
external magnetic field provide an example of a candidate detector technology that might 



be used. Consideration of detector technologies and their performance is under study 



and is outside the scope of the present paper. We will take M = 5 kt as a plausible, but 
aggressive, choice for the detector mass. 

Consider first an intensity-upgraded ENuF, namely a 20 GeV storage ring in which there 
are 10 20 muon decays per year in the beam forming straight section. The r appearance rates 
from v e — > v T and z/ M — > v T oscillations are shown as a function of sin 2 2# 13 and 8m\ 2 in Fig- 1§ 
for a 5 kt detector at L — 2800 km. In contrast to the — > v T background rate, which 
is independent of sin 2 26?i3, the v e — > v T signal rate increases linearly with sin 2 26*13. Note 
that for the LAM scenario with sin 2 26*13 = 0.04 and 5m\ 2 = 0.0035 eV 2 , there are about 10 
signal events per 10 20 /i + decays, and 300 — > v T background events. Thus it is desirable 
that the r sign mis-determination be less than of order 1 in 100. Whether this requirement 
can be met by placing hybrid emulsion or liquid Argon detectors in a magnetic field, or by 
the development of new v T detector technology, remains to be seen. 

Next consider the dependence of the sin 2 2#i 3 reach (for detecting 10 v e — > v T events) on 
L and the storage ring energy. Fixing 8m\ 2 = 0.0035 eV 2 , the sin 2 2# 13 reach is shown in 
Fig. ^ to improve with energy, and to be almost independent of L over the range considered 
except at the highest energies and longest baselines, for which the reach is degraded. For 
L ~ 3000 km, an energy upgrade from 20 GeV to 50 GeV would improve the reach by about 
a factor of 5. The energy dependence of the muon intensity required to meet our v e —>■ v T 
discovery goal is summarized in Fig. |j. We conclude that a neutrino factory consisting of 
a 20 GeV storage ring in which there are 10 20 muon decays per year in the beam forming 
straight section would enable the first observation of v e — > v T oscillations in LAM, SAM, 
and LOW-type scenarios with sin 2 2^13 > 0.01 provided a 5 kt detector with good r signal 
efficiency and charge-sign determination is practical. 



VI. MEASURING THE CP NON-CONSERVING PHASE 

In the event that the solar solution is LAM, CP non-conserving effects may be large 
enough to allow a measurement of the MNS phase 5 at a high-intensity neutrino fac- 
tory [|^,^,^]. The total rate of appearance events very strongly depends on the value 



of sin 2 26*13 - For example, Figs. |8] and |9| show the event rates versus L for the LAM solution 
in Table | with 5 = and sin 2 26 13 = 0.04, 0.004, and 0. For L less than about 5000 km 
the event rates for the BIMAX solution are about 25% higher than the sin 2 26*13 — curve. 
Although the rates decrease significantly with decreasing sin 2 26*13 , even for sin 2 2^13 = 
there is a residual signal from the sub-leading oscillation in the LAM scenario which may 
be detectable. 

Figures [Tl], and [TT] show predictions for the CP dependent ratio R efJL versus the baseline 
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L for a 50 kt detector in the LAM scenario with 8m\ 2 = 3.5 x 10 3 eV 2 . The error bars are 
representative statistical uncertainties. These figures present the results of calculations with 
decays/year and sin 2 26 n values of (10 20 , 0.04), (10 21 , 0.04), and (10 21 , 0.004), respectively. 
Results for phases 8 = 0°, 8 = 90°, and 8 = —90° are shown in each case, for both positive 
and negative values of 8m\ 2 . For these values of sin 2 2# 13 the event rates show a strong 
dependence on the sign of 8m\ 2 (due to different matter effects for neutrinos and antineutri- 
nos), and a smaller dependence on the CP- violating phase. Figure ^ shows the results for 
10 21 muons and sin 2 26^ = 0.04 in finer detail. Note that at small distances R efJL is not unity 
for (5 = (even though matter effects are small) because the and CC cross sections are 
different. We note that the CP-violating effect is largest in the range L ~ 2000-3000 km, 
vanishes for L ~ 7000 km, and is nonzero but with large uncertainties for L > 7000 km. 

Similar calculations show that for the SAM and LOW model parameters R efl is essentially 
independent of 8, verifying the conclusions of Sec. Ill that CP violation is negligible in these 
scenarios. The effect of matter in the SAM and LOW scenarios, which depends on the sign 
of 8m\ 2 and the size of sin 2 29 13, is similar to the LAM case. 

A nonzero sin 2 2# 13 is needed both to determine the sign of 5m 2 2 via matter effects, and 
to have observable CP violation from the sub-leading scale; however, whether matter or CP 
violation gives the largest effect depends on the size of sin 2 29i3. This is illustrated in Fig. [13|, 
which shows the ratio of wrong sign muon events versus sin 2 2#i3 for our representative 
LAM solar solution with 10 21 muons/year. For larger values of sin 2 2#i3, say above about 
0.001, the sign of 8m\ 2 (through matter effects) has the largest effect on the ratio, while for 
sin 2 29 13 < 0.001 the value of 8 (which largely determines the amount of CP violation) has 
the largest effect. 

Now consider the neutrino factory energy and intensity needed to begin to probe the CP 
phase 8 in the LAM scenario. We will define the sin 2 29\3 reach as that value of sin 2 2#i3 
that (with a 50 kt detector and two years of data taking) will enable a 3a discrimination 
between (a) 8 = and 8 = tt/2 and (b) 8 = and 8 = —ir/2. The measurement will be 
based on a comparison of wrong-sign muon rates when respectively positive and negative 
muons are alternately stored in the ring. The sin 2 2#i3 reach when there are 10 21 decays is 



shown for the 3a discrimination between 8 = and ±7r/2 in Fig. [14] as a function of baseline 
and stored muon energy. The optimum baseline is about 3000 km, for which the sin 2 2#i3 
reach is a little better (worse) than 0.01 for the 8 = it/2 (8 = — vr/2) discrimination, and 
is almost independent of muon energy over the range considered (Fig. f|). Thus, a high 
intensity 20 GeV neutrino factory providing O(10 21 ) muon decays might begin to probe 
CP violation in the lepton sector if the LAM scenario is the correct description of neutrino 
oscillations, and a 50 kt detector with good signal efficiency is practical. This conclusion is 



consistent with results presented in Ref. |3(| in which global fits to the measured oscillation 
distributions have been studied to determine the precision with which 8 and sin 2 29 13 can 
be simultaneously measured at a neutrino factory. 



VII. SUMMARY 

We briefly summarize the results of our study of the physics goals of an entry-level 
neutrino factory as follows: 
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(i) An entry-level machine would make a first observation of v e — > z/ M oscillations, measure 
the corresponding amplitude sin 2 2# 13 , and determine the sign of 5mJ 2 . 

(ii) The sin 2 2#i 3 reach for the first observation of v e — > oscillations and the measurement 
of the sign of £m 2 2 is insensitive to the solar neutrino oscillation solution (LAM, SAM, 
or LOW) within a 3-neutrino framework. 

(iii) A 20 GeV neutrino factory providing 10 19 muon decays per year would enable our 
entry- level physics goals to be met provided sin 2 28 ' 13 > 0.01 and a detector with good 
muon charge-sign determination and a mass of 50 kt is practical. The required beam 
intensity might be a factor of 2-3 higher or lower depending on where within the 
SuperK range the 5m| 2 parameter sits. The event rates also depend on the muon 
energy detection threshold (E min = 2-5 GeV) and we win or lose a factor of 2 in rates 
depending on how low this threshold can be pushed. To determine the sign of <5m 32 a 
long baseline (L > 2000 km) must be chosen. 

(iv) A candidate for an intensity-upgraded neutrino factory would be a 20 GeV facility 
providing 10 20 muon decays per year. In addition to the precise determination of the 
oscillation parameters, the upgraded neutrino source would enable the first observation 
of v e — > v T oscillations provided sin 2 2#i 3 > 0.01 and a 5 kt detector with good r signal 
efficiency and charge-sign determination is practical. 

(v) With a high-intensity neutrino factory providing a few x 10 20 muon decays per year, 
the ratio of [i + / ' wrong-sign muon rates might enable detection of a maximal CP 
phase in the case of the LAM solar solution. If sin 2 2# 13 is vanishingly small, (for 
example the bimaximal mixing scenario) a first observation of v e — > oscillations at 
a high intensity neutrino factory might provide a direct measurement of oscillations 
driven by the sub-leading 5m 2 scale, although wrong-sign muon backgrounds might be 
problematic. 

In conclusion, the required number of muon decays per year to achieve the various 
physics goals of interest are summarized in Fig. ^. At a 20 GeV neutrino factory ~ 10 19 
muon decays are required for a thorough search for v e — > appearance, ~ 10 20 decays to 
search for v e — + v T oscillations, and ~ 10 21 decays to probe the sub-leading oscillation scale 
and detect CP violation effects in a three-neutrino LAM scenario. 
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FIG. 1. Two patterns of three- neutrino mass spectra that can explain the atmospheric and 
solar neutrino anomalies: (a) <5m§ 2 > 0, (b) 5m| 2 < 0. 
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FIG. 2. Predicted measured energy distributions (including the detector resolution function) 
for v e — * CC events tagged by wrong-sign muons when 20 GeV positive muons are stored in 
the neutrino factory and a 50 kt detector is at L = 2800 km. The distributions are shown for a 
variety of muon detection threshold energies -EWn- The oscillation parameters are those for the 
LAM scenario (Table I). 
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FIG. 5. The ratio R of v e — * Vn to v e — > Vn event rates at a 20 GeV neutrino factory for 
5 = and ±7r/2. The upper group of curves is for 5m^ 2 < 0, the lower group is for 5m^ 2 > 0, 
and the statistical errors correspond to 10 20 muon decays of each sign and a 50 kt detector. The 
oscillation parameters correspond to the LAM solar solution with l^m^l = 3.5 x 10~ 3 eV 2 and 
sin 2 20i3 = 0.04. 
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FIG. 9. The same as Fig. || for wrong sign \x events from a /x + beam. 



1000 



100 



0.001 



10 zl n 

sin 2 29 13 = 0.04 



8 = 0° 
8 = 90° 
8 = -90° 




2000 4000 6000 8000 10000 

L(km) 
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FIG. 13. The ratio of v e — > to v e — > CC events versus sin 2 2#i3, with L = 2900 km, 
E„ = 20 GeV, 10 21 muons, and a 50 kt detector. The other oscillation parameters are the same as 
the LAM scenario in Table Q, and results for both positive and negative 5m\ 2 are shown. Predictions 
for maximal CP phases 5 = 90° (dashed curves) and 5 = —90° (dotted) are compared with the 
CP-conserving case 5 = 0° (solid). The error bars show typical statistical uncertainties on the 
measurements. 
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sin 2 2i>, 3 yielding 3a diff. for 6 = +-90° cf 0° 




FIG. 14. Reach in sin 2 29is that yields a 3a discrimination between (a) 5 = and 7r/2 with 
bm\ 2 > 0, (b) 5 = and 7r/2 with bm\ 2 < 0, (c) 5 = and — 7r/2 with 5m\ 2 > 0, and (d) S = 
and — 7r/2 with ctm 2 ^ < 0. The discrimination is based on a comparison of wrong-sign muon CC 
event rates in a 50 kt detector when 10 21 positive and negative muons alternately decay in the 
neutrino factory. The reach is shown versus baseline for four storage ring energies. The oscillation 
parameters correspond to the LAM scenario. 
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